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Abstract 

We propose a search strategy for the light R° (glueballino) particle sug- 
gested by G.Farrar in connection with the light gluino scenario. The basic 
idea is to moderate and stop the R° particles and then observe their decay 
to almost monochromatic 7r°-s - at an appropriate time delay relative to a 
primary collision event, where a gluino jet, likely to fragment into the R°, was 
produced. This technique is optimized at colliders (pp, e + e~), and depends 
on qualitative features of the R° hadronic interactions which we discuss in 
detail. 
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1 Introduction 



Over the last few years a case for light gluino and photino has been made by G.Farrar 
07 @> Hi- Supersymmetry breaking is expected to generate masses in the 100 — 10 3 
GeV range for the superpartners of the known light fermions and bosons. (The upper 
range is a theoretical bias and is natural if SUSY and E.W. symmetry breaking are 
concurrent and related — the lower limits are experimental "bounds"). However 
in certain schemes the partners of the exactly massless vector bosons (-associated 
with unbroken SU(3)c x U(1)em gauge symmetric), namely the photino and gluino, 
obtain only small <0(GeV) masses via radiative (loop) correction. 

Since finding a light gluino g at existing, pre-LHC, accelerators would have such 
monumental consequences, we believe that even the remote chance that this scenario 
is realized in nature is worth pursuing. 

Light gluinos modify various aspects of perturbative QCD - such as the running 
of a s (Q 2 ), quarkonium and jet physics ||. It seems that these modifications cannot 
- at present - definitely rule out the light gluino hypothesis. 

The gluino cannot, because of color confinement, directly manifest as a free 
particle, but rather as a constituent of a color neutral hadron the "glueballino" 
gg. This new hadron, termed by Farrar, is the lightest hadron with negative 
R parity. It is therefore strong interaction stable but decays weakly via squark 
exchange into a photino + hadrons 

R°^^ + H° (H° = 7T°, 77°, 7T + 7T-, TT + TT-TT , etc) (1) 

{If rrifjo < rrij the R° would be stable altogether. This theoretically (even more) un- 
likely scenario is forbidden by the following consideration: after the R°R ^hadrons 
annihilation freeze out, a relic density of R°: ~ 10~ 18 remains. This implies 
> 10~ 10 . As we will argue below R° is likely to attach to heavy nuclei. The 

Tip 

(R° — A, Z) composites constitute new exotic "isotopes" - which have been excluded 
with very high precision |7]].} 

The R° decay lifetime scales as the squark mass to the fourth power and is even 
more sensitive to the R° mass (or rather R° — 7 mass difference). For R° masses in 
the range 

m R0 = 1.5 ± 0.2 GeV (2) 

and 

70 Gev < nig < 500 GeV (3) 
Farrar (and we) find a very wide range for R° lifetimes 

r R0 ~(3x KT n sec - 10" 4 sec) (4) 

The R° hadron is a most striking prediction of the light gluino scheme. How 
come such a particle has not been discovered yet? 

As correctly pointed out by Farrar [|TJ, ||, searches for SUSY particles looking 
for missing Ej- (transverse missing energy) signals are rather insensitive to R°-s. 
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Over most of the above lifetime range the R° interacts in the calorimeters and loses 
most of its energy there. The final decay photino would then carry a tiny (< few 
GeV) energy and would be indistinguishable from a neutrino with similar energy. 
Likewise the R° could not have been discovered in beam dump experiments looking 
for penetrating particles. Finally searches for gluino jets based on their q jet like 
angular distribution || may also be somewhat hampered by the fact that the leading 
hadron in these jets, namely gg = R°, can take most (60%-80%) of the jet's energy. 
Unlike leading gg glueball in gluon jets which decay to multi 7r-s, K-s hadronic, 
10 -24 sec, time scales, the R° are stable on jet evolution time scale. Gluino jets may 
therefore be less conspicouse than gluon or quark jets. 

In principle neutral beam experiments with decay paths at a distance L from the 
production point can be used to look for R° decays if L/c ~ Jrot r o. Also the R°-s 
in the beam could be looked for via their hadronic collisions by careful measurement 
of time of flight and of R°p elastic collisions kinematics which hopefully can separate 
the R° from the dominant neutron component ||], 0- 

We will focus on an alternative approach. The hadronic interactions of the R° 
can be utilized to moderate it. R° decays delayed by t(R°) from "relevant" primary 
interaction, (chosen on the basis of being likely to send an R° in the detector's di- 
rection) could then yield almost monochromatic 7r°-s, ?7 -s and other 7r + 7r~, 7r + 7r~7r° 
hadronic systems. 

In the following we will review the basic steps and the features of the putative 
R° hadrons which control them. 

2 R° production 

The sensitivity of R° searches depends on its production rates and on the R° inclusive 
spectra, both of which we will estimate in the following. 

The gg-gluon coupling is larger than that of gg-gluon by the ratio of the second 
Casimir operators for the octet and triplet representations of SU(3) C : 

9m = 9 2 ggq = [C 2 (8)/C 2 (3)] g\ qg = 2.25 g\ qg (5) 

This suggests that gluino jets and mini-jets and the R° particles, that these jets 
fragment into, will be more copious than charm jets and charmed particles. This 
is even more so when the R° mass is appreciably smaller than the D — D* mass: 
Id = \m D * + \m D = 1.97 GeV. [m#o < tud is assumed to always be the case]. It 
was found that at the tevatron collider 10% of all jets had D + D* ||. We expect 
at least as many R°-s containing jets. 

Unlike the case of mesons, whose complete production pattern can be pre- 
dicted from perturbative QCD, a substantial fraction of D-s may be produced 
"softly" with no separate cc jets. This is a-fortiori the case for R° particle. If 
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this soft component can be described by the phenomenological fit 

da 



we would have (for m^o ~ 1.5 GeV) a fairly large ratio 

a(R°)/a(D + D*) = m * e ~^*° ~ 16 (7) 

Unfortunately 10 — 10 3 more neutrons and neutral kaons (Kq specifically) are pro- 
duced and provide a strong background for interacting (and decaying) R°-s, respec- 
tively. The R° particles are produced with larger transverse momenta than the 
latter n-s, K°-s and the other light hadrons. This feature is particularly true for 
particles produced in the forward (p) and backward (p) fragmentation region and in 
the central rapidity plateau (in the pp cms). This in turn causes a broader angular 
distribution of the produced R° particles 

p T (R) p T (n) . Pt(K l ) 
®r ^ — rBT > — rv or — — — ~ 6(n) or Q{K L )\ 8 
p{R) p{n) p{K L ) 

In particular many of the neutrons and some K -s may be "leading" particles, in the 
fragmentation of the proton projectile in a fixed target experiment. These will than 
have both small pt-s and large laboratory energies and hence very small angles. 

The neutral beams in various fixed target experiments, obtained by strong colli- 
mation in almost exactly the original proton beam direction, will therefore be further 
strongly enriched in neutrons and Kq. It is thus particularly difficult to look for the 
decays or interactions of the relatively rare putative R° particles in such beams. 



3 The lifetime of FT 

The most important parameter for determining optimal search strategies for the R° 
particle is its lifetime r^o. Indeed as tro spans the six-seven decade range indicated 
in eq.(j^) its very evolution as it propagates in matter will drastically change. For 
t r° — 10~ 9 sec it is likely to decay within a meter or so from the collision point and 
manifest via a missing signal. If tro > 10~ 8 sec then it is likely to become non- 
relativistic and further loose energy by elastic i?°-nuclei collisions before decaying. 
Let us therefore briefly recall why these lifetimes are naturally expected (if m^o is 
indeed in the m^o = 1.5 ± 0.2 GeV range prescribed!). 

The gluino decay g° — > qq^ with qq = uu, or dd, proceeds via the squark exchange 
diagram. 

If the gluino was a physical particle with mass rrig we could use this Feynman 
diagram to compute directly the g decay rate. By comparing with the W~ exchange 
in /i — > z/^z/ge decay we have: 

T(g) = a Q cD{m-g) a em tm W \i 5 /m 9 -x5 r(m|/m?) 

T(/i) (awfc) 2 ^ mq ' 3 ^mj r{ml/m%) 1 
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where the above four factors represent coupling constants and propagator ratios, 
N c (q 2 + q 2 ) = 3.5/9 = 5/3 is a charge color factor summing over uu, dd states of 
various charges, and a phase space ratio. (In particular 



r[z) 



]z + 8z 3 - z 4 -12z 2 \nz (10) 



is relevant in the limit 3> mj, Tn^). 

If the gluino had a "large", m~ ^ > 1.2 GeV, bare mass and, like the charm quark 
in D + = cd, dominated the R° = gg mass then, just like in the case of D + decay, 
we could use, in a spectator approximation, m g ~ m R o to estimate t(R°). Eq.(^) 
with o:qcd (1-5 GeV) ~ 0.3, a em ~ 0.0073, aw = 0.034 and <C m R o ~ m g (so 
that (m|/m^) ~ 1) yields: 

, n , r / "Tlo \ 4 / 0.1 GeV \5 
r(p) ~ r i?° ~ 1.6 x 10" 6 sec ( — — 9 - - ) ( ? -) 

~ 2.5 x 10" 12 sec (m 9 -/100 GeV) 4 (11) 

where we used in the last step ~ m R o ~ 1.5 GeV. However, in the envisioned 
scenario m~ 0; <C m^o, and eq.(|ll]) overestimates the R° decay rate. Most of the 
R° mass is due to QCD dynamic effects: the chromoelectric (and magnetic) fields 
i.e. surrounding gluon cloud and some qq sea. The bare gluino therefore carries a 
fraction xm R o of the R° mass with a (differential) probability f(x) normalized to 

i 

f(x) dx = 1 (12) 



o 



In a "generalized spectator" approximation we assume that the valence gluon and 
other gluons and/or qq pairs "sail along", as the bare, locally coupled gluino, of 
"mass" xm R o, decays. This then yields 



HR )- 1 = 0.4 x 10 12 (sec)- 4 r 100 GeV l 4 . f f ^ ^ r \( y / x f\ dx (13) 

L m g J J J 
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The ratio 



y ee (14) 
m R o 

controls further phase-space suppression doe to finite photino mass via r [(y/x) 2 ] = 

2 

7 — 1 -y2 . It is convenient to parameterize the decreasing probability that the 
gluino carries most of the total mass by power-like vanishing of f(x) as x — > 1 

f{x)~(i-xn 1 +i) (15) 

Even for = the last factor in eq.(|l3|) above supplies sizable suppression: 
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A(0) = / /(*) * 5 dx = (16) 

J 7 (0) is smaller than ~ 0.01 for 7 > 3. Clearly I 7 (y) = /(1-x) 7 x 5 r [(y/x) 2 ] < 7 7 (0) 
could be arbitrary small (or zero!) once m' — > m^o i.e. y — > 1. Following Farrar we 
will adopt the range of lifetimes: 

Ti?0 = (10~ 10 - 10~ 7 ) (rn g /100 GeV) 4 sec (17) 

corresponding to 2 x 10 -5 </<2x 10~ 2 . 

In the light gluino scenario the lower experimental bounds on squark masses, 
which are based on missing energy searches, do not apply. The squark decays (in 
~ 10~ 26 sec) to a quark and gluino. The gluino, after forming the R°, deposits most 
of its energy. 

For t(R°) estimates we will therefore take squark masses in the range 

60 GeV < m~ g < 500 GeV (18) 

The lower limit is based on the modification of R = cr e + e -^ hlldrons /a e + e -^ f j j +^- 
measured at LEP II due to qq production. The upper bound is a theoretical bias. 
Substituting in eq.(|T7|) we find that the wide lifetime range 3 x 10~ n sec < t r o < 
10~ 4 sec is indeed theoretically expected. 

The pattern of photino couplings implies that the qq state in the g — > qq + 7 
decay is 

2_ 1 - 

q u uu + qadd = -uu dd (19) 

3 3 

or, if we include the s quark, 

2 _ 1 - 1 _ 
q u uu + qddd ~ -uu — -dd — -ss (20) 

This state overlaps more strongly with the 7r° state 

71-° = (uu -dd)/V2 

than with i]° 

r]° ~ (uu + dd — ss) / v2 

In the absence of any dynamical "form factor" , effects this would enhance the R° — > 
7r°7 over R° — > rf^ by at least a factor of nine, i.e. 

T(R° -> r°7) /T(r]° -»■ r°7) ~ 9 (21) 

However, if m^o — m 7 ~ 1.4 ± 0.3 is on the high site of its range, it would appear 
that higher multiplicity final states R° — > 7r + 7r~7, 7r + 7r _ 7r 7 would be preferred 
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over both R° — > 7r°, 77 , +7. A crude estimate for the probability of the two body 
decay modes can be inferred from 



-7T S 3% „ ,„ o 



Br(D + ^ K ^) ~ 7% ~ 5r(i? u ^ 77r u ) (22) 



We will assume that 



7T° 



Erfi? -f 7 J = 0.1 ±0.05 (23) 
rf 



4 Hadronic i?° interactions 

An R° hadron propagating in matter scatters from the nuclei. These interactions 
are elastic if the kinetic energy of the R° is below the binding energy of nucleons 

Tro < B.E. < 10 MeV (24) 

and quasi-elastic [i.e. nucleons knock out] if: 

10 MeV < T R o < 0.5 GeV (25) 

Finally at higher energies 

Tro > 0.5 GeV (26) 

we have, in addition to the (quasi-) elastic -R°iVscatterings, also genuinely inelastic 
scattering involving production of pions, kaons etc. Most of these interactions are 
soft - at least in the t channel - and cannot be reliably calculated via perturbative 
QCD. Yet a qualitative understanding of these interactions and, in particular, the 
degree to which these differ from pion, kaon and nucleon interactions is required, if 
we wish to estimate the evolution of energetic (say E R o ~ P R o ~ 20 GeV) R°-s in 
matter. 

At very high energies it is believed that meson and nucleon scattering are dom- 
inated by the exchange of (two or more) gluons. At lower (Laboratory) energies 
-^Lab < 10 GeV quark anti-quark annihilations and exchange play a dominant role. 
In pre-QCD nomenclature these are the "Pomeron" and mesonic "Regge trajecto- 
ries" exchanges respectively with corresponding cross section: 

o" P om~giue-ex ~ const ( or InW ), W -> oo (27) 

coiisti' 

^Regge-quark-ex — ^ , W —> OO (28) 

What are the corresponding contributions for the case of R°N scattering? The 
fact that gluino (and gluon) have stronger gluonic couplings than quarks (eq.(|D 
suggests that the gluonic exchange (pomeron) part is enhanced there. Indeed the 
contribution of the chromomagnetic-electric cloud energy 

I R0 = J d 3 r (E 2 + B 2 ) [In R° state] - j [In vac] ~ m R o - m§ (29) 
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makes up most of the R° mass. Irq is larger than the corresponding integral in a 
meson or even in a nucleon. This conforms to the enhanced contributions of gluonic 
cloud interaction to R° — N scattering. 

On the other hand we have no analog of the quark annihilation or exchange 
diagrams for R° — N scattering. [We could have only gluon exchange analogies of 
these diagrams and even these only for say R°R° scattering]. 

Most of M — N and NN scatterings for E^.^ < 10 GeV is not only merely 
accounted for by Reggeon exchange. There are further detailed hadronic models 
involving pion and other meson exchanges for the specific reaction mechanism in 
various inelastic channels. Since to lowest order there are no such pion (meson) 
exchanges contributing to R° — N scattering we expect the latter to be smaller, 
particularly for Ei ab (R°) < 10 GeV. Further arguments suggest that 

a(R°N) < cr(MN) (30) 

holds not only for Ei a b(R°) < 10 GeV, but also at higher energies. The issue of high 
energy M N — NN cross sections is not completely settled. It is possible that "soft 
pomeron" contributions due to multiperipheral pion exchange type models make up 
a sizable part of the total M — N cross section even at lab energies as high as 20-40 
GeV. 

A closely related issue is the ambiguity of which meson nucleon scattering cross 
section should we compare a(R°N) with. These MN cross sections depend on 
whether we use 

M = n = ud : a tot (n + ) ~ 30 mb, (31) 



M = K = us : -[a tot (K + p) + a tot (K-p)} ~ 20 mb, (32) 

or 

M = (f) = ss : (J t ot{4>) = 10 mb. (33) 

{The total <j)N cross section in eq.(|33|) is inferred from vector dominance in high 
energy <fi photo-production. Estimates of <7t{pp) or &t{wp) done in a similar way for 
p, w photoproduction yield a(p°N) ~ a(w°N) ~ cr^^N) as expected for the p and 
w which like the pion tt are made of light (uu ± dd) quarks.} 

The clear cut systematics: 

&nN > &KN > V^N (34) 

is somewhat puzzling from the point of view of the gluon exchange model as the 
QCD coupling are flavor independent and g\ ug = g\ dg = g? sg . One may try and argue 
that the slightly larger constituent mass of the s quark (~ 500 MeV) as compared 
with ~ 350 MeV for the u, d quarks, can cause - via reduced mass effects - the 
sizes and cross sections to differ. Even if correct, this explanation undermines our 
original argument for an enhanced gluon exchange contribution to R°N asymptotic 
cross section in the first place. The "constituent gluon" mass - if this is indeed a 
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useful concept - and a-fortiori that of the gluino, are presumably significantly higher 
[say: m g , nig = 0.7, 0.9 GeV], than those of the s quark. The naive modeling of 
non-relativistic g — g states would suggest that these larger constituent masses (and 
stronger attractive forces) make a much smaller state. The smaller length of the 
"color dipole" would then compensate for the ~ y^9/4 = 3/2 larger color charges. 
The systematics (|34|) could reflect mainly the weaker contributions of pion exchange 
in KN and (even more!) in cf)N scattering, suggesting reduced R°N cross sections 
as well. 

One final piece of evidence for small i?°-nucleon cross section comes from recent 
HERA experiment where a small "pomeron-nucleon" total cross section has been 
inferred. This could be relevant to our discussion if indeed the pomeron is related 
to two gluon states. 

In the ensuing discussion of R° evolution in matter we will still allow a wide 
range of variation for the inelastic R° — N cross section 

40 mb ~ 1.4a 1 (NN) > <7 J (i2°JV) > 0.7a (<j)N) ~ 7 mb For T R o > 0.5 GeV 

(35) 

and similarly for the elastic cross section 

15 mb > a El (R°N) > 3 mb For T R o > 0.5 GeV (36) 

We should keep in mind, however, that the lower values are theoretically more 
favorable. Another crucial parameter controlling the hadronic shower initiated by 
an energetic R° is the fraction of the initial R° momentum, xf{R°), carried on 
average by the final R° in an inelastic reaction: 

_R°(momentum p 3> m R o) + AT (at rest) — > N' + n pions + R°(pf) (37) 

In proton-proton collisions a final "leading" nucleon (i.e. proton or neutron) carries 
on average xf{p) — 1/2- This presumably reflects the tendency of a "diquark" in 
the initial nucleon to reappear as a constituent diquark in the final nucleon. 

A key to the evolution of R° in matter is that the initial gluino in R° (= gg) 
must be conserved in all collisions. This implies that in any reaction initiated by 
an R° particle there must be an outgoing R° particle containing the initial g. Since 
the R° is likely to contain also the "initial constituent gluon", this leading R° may 
carry even a larger fraction of the incoming momentum than an outgoing nucleon 

x F (R°) ~ 0.65 ±0.15 (38) 

Indeed a large xp is also consistent with a "hard fragmentation" of gluino jets into 
a "strongly leading" R° carrying a large fraction of the gluino jet momentum. Such 
a hard g — > R° fragmentation may be required in order to avoid the exclusion of the 
light gluino hypothesis by the analysis of 4 jet LEP events (see ref.f|). 
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At intermediate energies |eq.(|25|)l the R° "knocks out" a nucleon out of the nu- 
cleus. The fraction of kinetic energy retained by R° can be approximately obtained 
just from the two body elastic R°N collision kinematics: 

Pl (R°) +p 2 (N) ^p[(R ) +p' 2 (N) (39) 

The kinetic energy of the final nucleon T 2 = E' 2 — is given in terms of the 
invariant momentum transfer t = q 2 , q = {pi — p[) via 

T 2 = t/2m N (40) 

In terms of the R° — N c.m.s. momentum p* and scattering angle 8*, t is given by 

t = 2(p*) 2 (l-cos9*) (41) 

Expressing p* in terms of the initial kinetic energy of the incoming R°, we finally 
find (for (p*/m N ) 2 < 1) that: 

n = Tl - T[ = - 2m l mR ° (1 - cosPK (42) 

The retained energy fraction T[/Ti is thus: 

rpi irr 1 2m N m R o 

so that on average we have for m R o = 1.5 GeV: 

(n/Ti)sP-N = 1 - 2m " m n° (1 _ ^ ! _ .4 8 (1 _ ^7^) (44) 

(m N + m R o) 2 



If the c.m.s. angular distribution is isotropic then cos^* = 0. This indeed will be 
the case at the lower end of the range considered where the scattering is mainly in 
S wave. If we allow a small admixture of P waves then 



cos6»* = 2(sin<5i/sin<5o), (45) 
where 5 , 5i are the S, P wave phase shifts. Assuming 

< sin5i/sin5 < 0.2 (46) 

we will have 

(Tj/Ti) = 0.62 ±0.1 (47) 
The ratio of the corresponding momenta or velocities is then 



p' 1 /p 1 ~ ft/ fa ~ y/Ti/Tt ~ 0.79 ± 0.16. (4J 
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Finally we address the region of slow (Tro < 10 MeV) R° particles, which should 
be discussed in terms of _R°-nucleus (rather than i?°-nucleon) interactions. 

The evolution of R° particles in this stage critically depends on the existence of 
a bound (R° — A, Z) state of the R° and various (Fe, Ca, S, Si and O) nuclei and 
on the probability that such bound states in reactions like 

R° + (A 7)^1 [R°(A-l,Z)]+n \ 

R +iAZ) ^\mA-i,z-i)] +P j- (49) 

For the long wavelength slow R° the individual nucleons will effectively merge and we 
expect that its interaction with the nucleus can be described via a smooth, effective, 
potential: 

lfc(r) = V. -f\ = V (50) 
p(oo) p(oo) 

where p(r) is the nuclear density and p(oo), V« D (= Vo) refer to an ideal case of 
infinite uniform nuclear matter. 

To simplify our estimates we will approximate also finite nuclei by a "square-well" 
potential: 

V B p- {AiZ) {r)=V -Q{R{A,Z)-r) (51) 

with r the distance of R° from the center of the nucleus and R(A, Z) the nuclear 
radius: 

R(A, Z) ~ 1.2 A 1/3 Fermi (52) 
If the momentum of the incoming R° 



p=^2m R oT R o (53) 

is sufficiently low so that 

p R(A, Z) < 1 - 2, (54) 

we may assume that the R° nuclear scattering is dominated by S wave. The last 
condition is satisfied for the nuclei considered with 56 > A > 16 only if T R o < 10 
MeV, which is indeed assumed here. Also eq.(^4|) justifies the smearing of individual 
potentials of R° interaction with the various nucleons. 
The S wave phase shift is given by 

S = — p R + arctani = ■ tan(\/p 2 + 2mVo ■ R)\ (55) 

VP + 2m Vo J 

The single parameter Vq = V R o summarizes our ignorance of i?°-nuclear interactions. 
We make the reasonable assumption that 

Vbp<\Va\~ 50 MeV (56) 

where V\ is the attractive potential seen by a A particle in a hypernucleus. Indeed 
the R°, just like the A particle, is not affected by the Pauli principle and can migrate 
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towards the inner shells. However, the R° — N interactions - lacking meson (and in 
particular pion) exchanges - are expected to be weaker. 

Indeed ir exchange diagrams contribute via 3 body interactions, to the A-nuclear 
binding. We have no analog of this for the R° case. As we will argue below (towards 
the end of this section), the companion of the R° particle is likely to be heavy. 

In particular 

m R + - m R o > (57) 

as compared with 

m E - m A o ~ 70 MeV (58) 

suppressing the three body putative diagrams and V R o. 

The following values of V R o, consistent with the bound |Vro| < 50 MeV, lead 
to qualitatively different behaviors of i?°-nuclear scattering in the T R o < 10 MeV 
range: 

(i) V R o is positive or negative, but very small: 

|Vro| <Tro ~ 10 MeV (59) 

(ii) V R o is positive (repulsive) and 

V R o ~ Tup ~ 10 MeV (60) 

(iii) * V R o is negative (attractive) and of magnitude 

2-4MeV< |Vro| < 10 MeV (61) 

(iv) * V R o is negative and 

50 MeV > \V R o\ > 10 MeV (62) 

(v) V R o is positive and 

50 MeV > |Vro| < 10 MeV (63) 
We will next discuss each of these cases. 

(i) Expanding the expression ([55]) for 5o(p) in the small quantities uiVqR 2 and 
mVo/p 2 we find 

^-^(tanpfl--^) (64) 
p 2 cos 2 pR 

Thus for either sign of Vq we obtain fairly small R° nuclear cross section 

a\R«- Nucleus) - ^ = const (65) 
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(ii) In this case S can be appreciable 

aW(R° - N) ~ % ~ 7r(i? AiZ ) 2 , ^ = const (66) 

(iii) * The starred cases (hi)* (and (iv)*) are distinguished by two important features. 

First we are guaranteed to have one (or more) bound states. Indeed the 
condition for having at least one S wave bound state: 



/2mV R(A, Z) > tt/2 (67) 
amounts to 

Vro > 4 - 2 MeV (68) 
for nuclei in the oxygen-iron range i.e. for 

16 < A < 56 (69) 

Second we will argue later that cases (iii) or (iv) are the most likely in any 
event. 

For case (iii) we expect enhanced S wave scattering cross section, i.e. 

However the small binding energy (|e| < |Vo|) may kinematically disallow (or 
suppress) the actual R° capture reaction fl4"9] ) since there may not be sufficient 
energy for knocking a nucleon out. The latter requires in particular that 

T R o + \e\ > B.E. of nucleon ~ 8 MeV (71) 

(iv)* In this case we could have several S wave (and perhaps even P wave) bound 
states. However the most important point is that eq.(|7T|) is now satisfied, 
sometimes with a wide "margin" . 

The R° capture - reaction fl49D , leading to the formation of the R° — (A — 
1) bound state, becomes then exothermic and the capture cross section is 
enhanced by the ratio of outgoing nucleon momentum and the incoming R° 
momentum: 

p N (final) ^ j T R o + jej — B.E. m n > ^ 
Pro V T R o m R o ~ 

The total scattering cross section in this case af^ is larger than af^t ■ The last 
case corresponds to a flipped (repulsive) potential: V^ v ' = -7™. Hence (see 
sec.(v) next) 

4? > a? ^ nRf AZ) (73) 
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The kinematic enhancement of eg. fl72) suggests that the capture component 
in the total cross section 

= <#° + (74) 

is appreciable. 

(v) In this case the nucleus acts like a hard sphere. Furthermore pR(A, Z) > 1, 
and we may use the classical cross section 

aW = ttJ^,. (75) 

For elastic i?°-nucleus collision the energy fraction carried by the outgoing R° is 
given by eq.((43|), but with m N — > m(A, Z) ~ Am^r. Thus the average kinetic energy 
fraction retained 

2uA 

(Tl/TOHO-Nucleus * 1 ~ ^2 + 2 ^ + ^2 ' U " C0S ^) ( 76 ) 

with 

m = — — = 1.6 ±0.3 (77) 

is significantly closer to one that in the case of elastic i?°-nucleon (a feature well 
known from neutron moderation). Eq . (|76|) implies a corresponding ratio for outgoing 
and incoming R° velocities: 

uA 



<#/A>j*>_ Nucleus = 1 - A2 + 2uA + u2 ■ (1 - cos^) (78) 

We next proceed to estimates of Vq - making first an argument for negative 
i.e. attractive Vo- Approximating R° ~ g+gluons, i.e. neglecting additional qq 
constituents, there are no Pauli exclusion repulsive effects in R°N scattering. In fact 
qq annihilations and quark exchanges i.e. ordinary mesonic contributions to R° — N 
scattering, repulsive or attractive, are altogether absent. The only interactions are 
gluonic exchanges i.e. the QCD analogue of "Van-der-Waals type" interactions. 
Color confinement implies a finite lowest mass rn^g" > in the gluon exchange 
channel. This will modify the r -6 — r~ 7 potentials derived in the perturbative two 
massless gluon exchange approximation by an exponential factor 

e- m V (79) 

We conjecture however that the generic property of the 27, 2g perturbative exchange 
potentials, namely their almost universal attractive character, will survive in the 
full-fledged QCD treatment of R°N interactions. 

In the perturbative approximation the "Casimir- Polder" interaction between 
two hadrons A and B is proportional to the following combination of electric and 
magnetic polarizabilities of A and B II: 



C 



AB 



a W a ( B ) + a W J B ) _ (JA) J B ) + JA) Jfl)s 



JO) 
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For A and B hadrons which are ground states in the respective channels, second 
order perturbation in an E(B) background fields, say ~ a^E 2 , imply negative 
a l E (or a l M ). Hence Cab is guaranteed to be negative once at least one of the 
following holds: 

(i) Chromoelectric effects dominate: ^> a^j 

(ii) = affi and cuff = 4f, or at least 

AA) I (A) _ (B) i (B) fo, \ 

a E l a M - a E l a M \° L ) 



Hadronic constituents are generally more relativistic than valence atomic electrons 
so that chromo-magnetic effects are more important. Also three and more gluon 
exchange are more important relative to two photon exchange. Thus the attractive 
character of the R° — N gluonic interactions is not obvious. We proved [[10] via 
rigorous "QCD Inequalities" techniques that the low energy interaction between 
ground state pseudoscalars 

M ab = "g a7 5%", M' cd = u q c y 5 q d " (82) 

is always attractive if 

m a = m c and = rrid (83) 

but all flavors a, b, c, d are different. The latter condition does indeed exclude any 
q exchange or qq annihilation so that this MM' attraction refers specifically to the 
gluonic interactions of interest. 

We believe that the same result holds for the interaction between any two hadrons 
of similar color structure - reminiscent indeed of condition (|8l|) for the perturbative 
case. 

The gg = R° and N = uud nucleon states are clearly different. There is, however, 
no obvious reason why the sign of their mutual gluonic interactions be reversed 
relative to the gluonic interactions for iViV or R°R°. We will therefore assume a 
negative Vo- Unfortunately the magnitude |Vo| is not known. We have suggested 
above (eq.([5§)) that \V$ R) \ < |F (A) | ^ 40 - 50 MeV. We will argue next for the 
following likely lower bound: 

\V (BP) \ > 10 MeV (84) 

We express the fact that most of the R^ mass originates from QCD dynamics in 
the following suggestive manner. 

1 GeV < m R o - mf ~ / d 3 f (E 2 C + B 2 C ) (85) 

"over BP" 

In the presence of nuclear matter the last expression will be modified. This modifi- 
cation can be parameterized by the chromo-dielectric (magnetic) constants shifting 
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away from unity. Assuming that / E 2 dominates, we have for the rest energy shift, 
i.e. for Vo, the following expression: 

V = (~ - 1) / E 2 d 2 r = (i - 1) (m R0 - m°) 

~ -(e- l){m R o -mjj) < -(e - 1) GeV (86) 
where we used eq.(|35"D and in the last step expanded in 

e - 1< 1 (87) 
We note in passing that Vo < amounts to the natural assumption that 

e - 1 > (88) 
To estimate e — 1 we use the analogue of the familiar QED expression 

(e -l) = £„„«»> (89) 

with the (chromo-) electric polarizability of the nucleon and 

n = (1.2 Fermi)" 3 (90) 

the nucleon number density. Ordinary electric polarizability can be written as 

a E ^ r'l (91) 

where in the last expression ro is an effective size |eq. (|9~T|) also represents the classic 
polarizability of a conducting sphere of radius ro]. Collecting all the above results 
we find 

N,Eff n 

i r) GeV (92) 

.bermi ' 

Thus | Vo| > 10 MeV once the mild requirement 

r N,Eff ^ 1? Fermi ^ g3 ^ 

is satisfied. 

Before concluding this section we would like to reiterate the feature which is the 
most important distinction between the i?°-Nucleon scattering and the scattering of 
any other stable hadron on a nucleon. 

All the hadrons with lifetimes > 10 -10 sec (with the exception of A ) come in 
isospin multiplets: 



7T 
7T° 

7T~ 



'K+\ (K°\ fp" 
.K° \K- U, 




/ 5+ \ 
V 5- / 
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Also the A is associated in SU(3) flavor, with the S multiplet from which it is split 
by only 70 MeV (see eq.®). 

In principle we could have in addition to R^ = gg other more complex hadrons 
containing a gluino which could be charged like 

R + = g{ud) 8 (94) 

where as indicated the ud are in a color octet combination. However the R + is 
definitely not an I-spin partner of R°, but rather the lowest member of a completely 
different family of particles. These particles are the SUSY partners of the "Meikton" 
or "Hemaphrodite" hadronic states, which have been speculated by various authors, 
namely g(ud) 8 , just as R° = gg is the SUSY partner of the glueball state: gg. 

To date we have no clear-cut evidence for either glue-ball or "Meikton- 
Hemaphrodite" states. Theoretical prejudice suggests that the glue-balls are sig- 
nificantly lighter: 

m (gg)<m°[g(ud) 8 ] (95) 
which in turn suggests that also 

m R o = m°(gg) < m° [g(ud) 8 \ = m R + (96) 

Fortunately we have also direct "experimental" evidence that the R + should be 

heavier than the putative light R° by at least m^. (see eq. (|57D ). Indeed if the 
reverse, namely 

m R + < m R o + 7r + , (97) 

is true, then we would have a new, (almost) stable, charged hadron of mass m R ± < 
1.6 ±0.2 GeV. 

Its production cross section - particularly at large px - should considerably 
exceed these of the deutron d or antideutron d. Such particles should have been 
detected for lifetimes in the range considered 

3 x 10~ 9 < t r+ ~ tro < 10~ 4 (98) 

Thus the R° is really conserved in hadronic collisions. In any reaction initiated by 
R° the final state will also include an R°. In particular, any R + produced would 
decay via 

R + -> R\ + (99) 

on hadronic time scales. 

One final comment concerns the small nuclear effects on R° nuclear scattering at 
high and intermediate energies. These "shadowing" effects tend to reduce the cross 
section of any projectile x on a nucleus from the naive, weak coupling value 

a {0 \x,A) = A cr(x,p) (100) 
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and correspondingly prolong the interaction mean free path from 



/ (0) (x, A) = , = ——. (101) 

na{x } p) pAva{x,p) 

with Av = Avogadro number = 6 x 10 23 and p the material density so that n = pAv 
is a nucleons number density. The shadowing effects are appreciable when 

a°(x, A) = A a(x,p) ~ tcR 2 {A,Z) = ttA 2/3 (1.2) 2 (Fermi) 2 (102) 



or 



a{x,p) > A- 1/3 ■ 45 mb ~ (18 - 12) mb (103) 



for A = 16 - 56. 

Thus we will use for iron longer interaction path than those implied by eq.(^), 

my- 

lmfp{R° -Iron) = (10 - 30) CHI (104) 

with the lower (higher) value corresponds to &^ R o_ N ) = 40(7) mb with substantial 
~ 50% (no) shadow corrections. 



5 The evolution of energetic i? u particles in mat- 
ter 

Energetic, initial (-R ) > 30 — 40 GeV, R° particles evolve as they propagate in 
matter in three main stages corresponding - in reverse order - to the energy ranges 
of eqs.(p4l),(p5D and ( p6|) above. We will next try to follow the evolution in these 
stages - utilizing the estimates of the previous section. 

Stage (a): Fast R°-s are conserved through their hadronic cascade evolution - if 
we neglect the presumably small probability to generate further R°R° pairs in R°N 
collisions with c.m.s. energies Wron — (2-E/jomjv) 1 / 2 < 10 GeV. This drastically 
simplifies their evolution as compared with that of other long lived neutral particles 
such as K°-s, which are created in secondary hadronic collisions, and neutrons, 
which, even neglecting secondary iViV pair creation, can be knocked out of nuclei. 
All these processes not-withstanding the hadronic showers generated by qq jets (and 
also cc, bb jets) essentially quench at Z que ~ 10 hadronic interaction lengths inside 
the calorimeter or shielding materials. 

We would like to argue that as they evolve through the energy range 40 GeV > 
Tro > 0.5 GeV, the R°-s will typically get further away from the initial interaction 
point than l qne , i.e. further than typical hadronic showers. Several effects contribute 
to this: 

(i) Charged particles continuously loose energy via ionization and stop at well de- 
fined ranges - even if they have no hadronic interactions. These losses increase 
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at low (sub GeV) energies: thus the ranges for 2 GeV p, K + and tt + in iron 



arc: 



R(p) = 16 cm, R(K + ) = 30 cm, R(tt + ) = 60 cm (105) 

These ionization losses can indirectly contribute also to the slowing down of 
K°-s and neutrons. At T# ~ 1/2 — 2 GeV hadron's kinetic energies, we have 
appreciable cross sections for charge exchange reactions 

K°(K°) +N -> K + (K-) +N' (106) 

n + iV -> p + N' (107) 

with iV referring to a nucleon in the nucleus. Also in inelastic collisions at Tjj = 
2 — 40 GeV the incoming K° or neutron often "fragments" into a leading K ± 
or proton. Thus roughly during half of the evolution of the hadronic showers 
the nucleonic or "kaonic" component manifests as the charged component of 
the isospin multiplets which slow down via ionization. 

As emphasized towards the end of the previous section there are no analog 
R° — ^stable R ± conversion and there is no - even indirect - ionization losses 
during the R° evolution Jl2| . 



(ii) The R° particles are likely to maintain higher fraction xf{R°) of their initial 

momenta than nucleons or kaons. 

(iii) The i?°-nuclear inelastic cross sections in general and particularly in this 
40 GeV > T R o > 1/2 GeV range are likely to be smaller than N — N and 
even KN cross sections. 



To simplify the treatment of R° evolution in stage (a) we will neglect the slow 
buildup of transverse momenta - namely momenta transverse to the initial entry 
direction of the R° particle into the material. (After n collisions we have 

\pP\ 2 ~nA 2 (108) 

with A — 1/2 GeV typical soft momentum transfer.) This leads then to a one- 
dimensional evolution in momentum space [or more conveniently in rapidity space 
y n = — In x n with x n the fraction of the initial momentum carried by the R° after n 
collisions] and in coordinate space. The coordinate is the distance r (measured along 
the momentum of initial R°) from the entry point to the material, to the R° location. 
For energy independent cross sections (particularly appropriate for R°N scattering 
which are devoid of the varying Reggeon contribution) these y and r evolution 
decouple. In a uniform medium, assumed for simplicity, the mean interaction length 
for R^ is then a constant, 1$. Let gM(r) be the r distribution after n collisions. By 
definition 

gQ-\r) ~ e~ r/l ° (109) 



18 



We can readily derive the recursive relation 



g {n+l \r) ~ J dr' g {n \r') e - {r - r ' )/l ° (110) 
o 

and its solution 

g( n \r) ~ e~ r/l ° r"" 1 (111) 
The average distance after n steps is then 



r W = nZ = (112) 

as expected. 

The rapidity space evolution will be completely analogous if we assume that after 
a single collision we have rapidity distribution f^\y) = e~^ y (corresponding to x 13 
in x space) so that y^> = 1/(3. 

Let us denote by A' a ' the energy-momentum average reduction factor in stage 

(a) 

A« = x { F R0) (113) 

According to eq.(|3"8"|) we take A^ = 0.65 ± 0.15. Starting with an initial energy 
Eq = T (= 40 ± 10 GeV), we need then altogether N^ a > collisions in stage (a) so as 
to degrade the kinetic energy to the required Tf(= 1/2 GeV) where 

N (a) = ln(^o/T/) Q ± 20 w ln(L6 ± Q 3) ^ 7 _ 18) ( 114 ) 

ln(l/A) 

[We will consistently use () brackets following equations for the actual numerical 
estimates used in applying the equations.] The total penetration distance traveled 
by the particle in stage (a) is then on average 

L (a) ^ N (a) j(a) = ^ ± g^(a) = (3 ± ^ eqmvalent ^ ( U5 ) 



where we used eqs. (|114| ) and ( |104j ). The very large spread in ( |115| ) represents the 
accumulated uncertainties in a I (R°N) or l[^ t and in the elasticities X^ a ) and finally 
in the initial R° energy E {0) = 40 ± 10 GeV. We believe that the values = 0.7 
meter i.e. or = 5.3 meter i.e. are very unlikely. In particular to achieve L m i n 
we need to assume a^?]^ 3 * 10 ~ 40 mb (and xp(R°) ~ 1/2) which are too high (low) 
respectively and for L max we need a 1 ~ 7 mb (and xf(R°) ~ 0.8) which appear to 
be too low (high) respectively. We will adopt 

L (a) = 3± 1.5 meters (116) 

The main conclusion which we would like to emphasize is that the putative R° 
particle travels beyond the full extent of (e.m.+hadronic) showers due to normal 
(non R°) hadronic jets. Indeed the extent of the normal shower is affected also by 
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ionization losses of the leading particles to which the R°-s are completely immune. 
Throughout most of the stage (a) the R°-s move with (3 R o ~ 1, hence its total 
(Lab-time) duration is 

t (a) ^ L (a), 1Q -8/ 1 + 3 n / n7) 

' v -0.4 ; v ; 

Because of time dilation effects the required proper time (and hence the lower bound 
on r R o required if indeed the R° particle is to complete this stage) is lower than the 
above t^ a \ Specifically we find that 



i=0 C c 



A" 



~ (2-5)- ~ (0.6-5) x 1(T 9 sec (118) 

c 

where in the last equation, 7$, the proper time lapse between the — ith and 
N a —i + 1th collision, is shortened relative to Iq/c by I/7W ~ (A a ) 1 and we neglected 
(^a) Na = ^ — Thus the modest requirement 

r^o > 5 x 10~ 9 sec (119) 

suffices to ensure that the initial R° particles survive stage (a). 
Stage (b): Intermediate energies. 

Using eq. ([|7|) for the average fraction of kinetic energy A^(= 0.62±0.1) retained 
in quasi-elastic _R°-nucleus collisions in stage (b) we need on average 



JV« = - 



r (b) /T (b) 
initial/ final 



/,,rj0 3) (120) 



ln(l/AW) v In 1.6 ' 

collisions to complete this stage. We used in ( |120| ) T^ tial = 0.5 GeV, T^ al = 0.01 
GeV for the kinetic energies characterizing the entry to and exit from "stage (b)". 
The kinematics of the elastic i?°-nucleon collision implies that for Jfr(-R°) —constant 
(i.e. constant differential cross section in the cms frame) we have in the lab frame 
for u = m R o/rriN |l| : 



At l + w 2 cos2# 1 + (2.5 ±0.6) cos 26 

f u {cos6) = —\ R0 - {1 _ u 2 sin 2 9 y /2 ^ (i_( 2 .5±0.6)sin 2 ^)V2^ ( 121 > 

Let the initial R° direction be no, and h k be the direction after k successive scat- 
terings. In f u (cos9) the cos# argument refers to h k ■ h k ^\. The distribution after k 
steps will be 

/ (fc) (cos0) = f (k) (h k -h ) = J dh 1 dh 2 ...dh k 

xf u (n ■ 77,1) f u {h x -h 2 ) ... f u (h k -i ■ h k ) (122) 
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Using the completeness relation 

j-J dh' P l (n.n')P l ,(n'.m) = ^^5(l,l')P l (n-m) (123) 

we find that if 

oo 

f u (cos9) = E /,(«) (21 + 1) P,(cos0) (124) 
z=o 

then 

oo 

/<*>(cos0)=£ (2Z + l)fl(oosfl) (125) 

We will be particularly interested in the net extra displacement of R° - namely the 
extra displacement along the initial R° direction (no) accumulated during stage (b): 

AL^ — l Yl cos(n k no). On average it will be 
fc=i 

AW = /oE (n k • no) = k E (A)" = T^ 1 " fi) ( 126 ) 

fc=l k=l 1 Jl 



with 



1 + v 2 (2r 2 — 1 "l 

h(u) = / d(cos9) f u (cos9) P 1 (cos9) = C / dxx—^-± — (127) 

J ^ |(1 — u z ) + iran ; 

(_l +u 2)-l/2 



[(l- M 2 )+w 2 a; 2 ] 1/2 

where 

a 

and C is a normalization constant 

„ , /• , l + w 2 (2a; 2 -l) , x 

C- 1 = / rfx 1 — ^ ^ 128 

{ [(l-« 2 )+w 2 :r 2 ] 1/2 

The transverse displacement is given by a "random walk" 

(A?) 2 = (/f) 2 E sin 2 (n^no)(< |i\^Z 2 ) (129) 

using the above with ~ Z b ^ ~ 10 — 30 cm, u — 1.6 ± 0.3, A^ 6 ) = 8 — 3 we 
estimated: 

AL (6) ~ 0.4 - 1.3 meters i.e. (130) 

A^ 6) ~ 0.5 ± 0.25 meters (131) 

Finally the R° motion during most of stage (b) non-relativistic. The velocity after 
the n-th collision is therefore: 



& = a,(Vaw) b ( 132 ) 
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The total time (~proper time) consumed in this stage is 

At ( b ) = Ar W = E At (6) = |_ 1.5 - 5 x 10- 8 sec) (133) 

with /3o — 0.8 the initial velocity corresponding to Tro = 0.5 GeV. 
Thus if 

tro > r {a) + Ar {b) ~ (2.5 +6 ) x 10~ 8 sec (134) 

—0.5 

the R° particle would, on average, live through the end of stage (b). Its distance 
from the entry point is then: 

L ab = L {a) + AL {b) ~ 4 ± 1.6 meters (135) 

with a transverse spread of 

A T = 0.5 ± 0.25 meters (136) 

Stage (c): Our discussion in sec. 4 suggests that once T R o < 8 MeV (or the 
R° velocity satisfies (3 < /3 = 0.1) the R° particles will be quickly captured into 
R° — (A, Z) bound states. This would be effectively "fix" the R° and prevent further 
diffusion away from the R° location it reaches by the end of stage (b) at the (L, At) 
values quoted above (eq. (|135|) , (|136|) ). 

Even if the R° particles stay unbound for the duration of its lifetime t r o , the extra 
diffusion distance V N( C H grows very slowly with t r q like ^In t r° ■ From eq.([78D 
above we find that the time intervals between successive i?°-nuclear collisions are: 

AtW = ^( 1 + i)"^( 1 + l)" < 137 ' 

Thus all V^ c ^ collisions of stage (c) will be completed after a time 

7V< C ) A A 

t® ~ t b0 = V At^Ul + -) iV<c)+1 - l}- = At^-e^ (138) 

1 A > u u 

where 1 is implicit. Specifically we find for tro = 10 -4 sec 

N {c) = — ln (^^) (- 80 for A = 56 = 15 cm u = m R0 /m N ~ 1.6) (139) 

where Iq was estimated via the geometric R° — (A, Z) cross section o = ttR 2 (A) = 
3.6 x 10" 25 cm 2 to be 

Iq ~ 15 cm in iron (140) 
For r R a = 10~ 6 sec we find instead of eq. 



N (c) = 40 (for A = 56, l = 15 cm, At = 5 x 10" 9 sec) (141) 
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The maximum total extra diffusion distance added in stage (c) is therefore 



|A# (c) | = v/iV^Zo =\ lA meterS ' TR ° = 10 "! ' 

I 0.9 meters, r R o = 10~ b sec 



(142) 



In passing we note that once R° (or a neutron) achieves a final thermal velocity (~ 2 
km/sec, or (3f = 0.6 x 10~ 5 ) one reverts again to ordinary diffusion with AR ~ y/r. 
However even after 100 elastic collisions with iron nuclei a R° particle with starting 
velocity /3 = 0.1 will on average go down only to Pf = 0.005. 

A key observation is that if the R° decays at any time during stage (c) we expect 
its final velocity to be 

Pf < 0.1 = Po (143) 
Furtuitisely this holds also if the R° is bound to the nucleus with binding energy 



|e| ~ 8 MeV. The Fermi momentum of the bound R° (~ y2\e\m R o ~ 0.16 GeV) 
corresponds to 

P f ~ 0.1 (144) 

Consequently the Lab energy of the 7r° from the decay R° — > 7r 7 will be rather 
monochromatic. Specifically 

E n o = (l + frcos*)^ = (l + ftooBfl ™*"™^"^ (145) 

satisfies (on average) \E w o(f3,cos9) - E$\/E$ <% = 0.05. 

We conclude this section by a short discussion of the likely location of this 
R° decay event, r (decay), measured relative to the intersection point. We will 
consider first the longitudinal displacement tl along the initial R° direction, ho- It 
is accumulated only in stages (a) and (b) and assuming t r o > 10~ 8 sec so that stage 

(b) was completed we have from eq. (|135|) 

r L ~ L( )+(6) = L ab ~ 0.4 ± 1.6 meters (146) 

The total displacement has also random diffusive "components" from stage (b) and 

(c) . Thus finally the total displacement is 



if) 1 ' 2 ~ r ~ V ( L ab) 2 + (Ag) 2 + (i? c ) 2 (~ 4.5 ± 2 meters) (147) 

And its average direction relative to no is 

h • r ~ r L /r ~ 0.8 (148) 

The above estimates referred to iron, and also neglected the freely traveled dis- 
tance to the calorimeters and hence r is underestimated. We note that lighter (con- 
crete rock+earth materials) may actually expedite the slowing down of R° particles 
in stage (c) because the fraction of energy loss is ~ -r. 
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6 The proposed i? u search strategie 

The search strategy proposed here address a wide range of putative R° parameters. 
Specifically we will assume r R a > 3 x 10~ 8 sec, cr R a N < cr NN for both elastic/inelastic 
cross sections and xp(R°) > xp{p) = 0.5 for the momentum fraction retained by the 
leading R° particle in inelastic collisions. The search is optimized if <j r on — ^&nn, 
xf(R°) ~ 0.7 and particularly when Br(7r°), the branching ratio for the R° — > 7r°7 
decay mode is large, say Br(n°) > 20%. 

Let us assume that a detector "D" of volume Vd has been installed underground 
at, say, the FNAL collider (or at e + e~ colliders) at a distance _R(=5±2 meters) from 
the intersection point and with R perpendicular to the beam direction. This detector 
searches R° — > 7r°7, rf^ and possibly also R° — > 7r + 7r~7, R° — > 7r + 7r~7r°7 decays. 
The signature distinguishing such decays from background events could be an almost 
monochromatic 7r° with E n o > 0.6 GeV. Alternatively final states including several 
pions (and possibly some nucleons knocked out when the decay R° — ► 7r + 7r~7r°7 
occurres inside the nucleus, once the (R° — A, Z) bound state formed) should have 
a relatively isotropic distribution with limited missing momenta p^ < m R o/4 ~ 0.4 
GeV. Our arguments for a longer migration length for R°-s as compared with the 
extent of hadronic showers suggest a reduction of the neutron background at the 
location of the detector if it as placed at a distance r given by eq.( |147| ) above. 

The collision of such neutrons with nuclei in our detector can produce ir°, 7r + 7r~, 
etc. final states and constitute a background to the R° decays. However only 
neutrons with kinetic energies T R o >GeV (T R o > 2 GeV) are likely to produce n° 
(7T + 7r~7r°) in R°N collisions. These background events will have achromatic tt° or 
apparent missing momenta - often pointing in the direction of f i.e. towards the 
decay vertex. 

One can further reduce the energetic neutrons background in our detector by 
excluding time "windows" At, say of magnitude At ~ 3 x 10 -8 sec following each 
collision event or, if the rep cycle of the collider is short (as in the upgraded FNAL 
Tevatron), the first At ~ 3 x 10~ 8 sec following relevant events with jets. [Only 1-2 
GeV neutrons from such events emitted in the direction of the jets are likely to make 
it to the required transverse distance from the colliding beam. In typical minimal 
bias events the neutron along with most other particles produced in the beam's jets 
are forward backward peaked]. The total distance traveled by relativistic hadrons in 
this At is 9 meters, and most of the energetic neutrons are likely to have interacted 
and slowed down by then. Yet the slowly diffusing R° may well have not decayed 
if t r° > 3 x 10~ 8 sec. This negative time correlation will also reduce background 
from energetic muons emerging from the primary collision which interact in the 
vicinity of our detector and produce secondary neutrons. Such muons, cosmic ray 
muons, and charged particles background in general, can be further reduced by 
an anti-coincidence network surrounding our detector at a distance of ~1 meter 
from its boundaries. Also the very strong, slow or even thermal, ambient neutron 
background - causing many low energy depositions, via nuclear captures - can be 
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reduced by using a proper neutron absorber shielding. 

Finally we note that having over most of the distance f between the intersection 
point and the detector iron magnetized in a direction perpendicular to r could also 
drastically reduce the neutron background by having the intermediate proton curve 
in the magnetic field [T2] . 



All this notwithstanding we cannot hope to conclusively "pin down" R° particles 
by looking at their decays only. The basic reason is that, even if the detector is placed 
at an optimal distance r (see eq. ( |147| )) (which is a priori known only within ±60%!), 
and the R° particles are isotropically distributed, only a small fraction 

/ = VD Br[ f^ ^ * [V D /m 3 ] • (3 x 10- 5 - 2 x 10- 3 ) (149) 

3 r 

of the R°-s will decay via the 7r°7 mode, inside the detector. [In the above estimate 
we used r = 5 ± 3 meters, Br = 0.1 ± 0.05]. 

To pick out the genuine R° decays we correlate events in the detector with 
primary pp collision both temporally and directionally. Specifically we will focus on 
primary pp collisions with g jets that could give rise to a sufficiently high energy R° 
that in turn could make it to our detector. 

Thus we will consider those "special" pp events with one (or two) jet(s) of trans- 
verse energy Ex > E^\~ 10 GeV), and with no energetic lepton(s) [to avoid b,c 
jets] in the jets(s). If a light gluino exists these would be gluino jet(s) with an ap- 
preciable ~20% probability. Beam crossing and one pp collision on average happens 
every 

S tg = 3 x 10~ 6 (10- 7 ) sec 

in the old (upgraded) versions. However only one in Njj ~ 10 6 collisions is "special" 
in the above sense. (N H is the ratio of total pp cross section to a perturbative 
estimate for the cross section of the gluino-jet production). 

Hence the "special" collisions will be spaced, on average, by Ai# = NnStn — 
3 — (0.1) sec. If the AE > 0.6 GeV energy deposition in our detector is due to a 
decay of an R° from a gluino jet in a particular "special" event then we expect that: 

(i) The decay should happen on average at time ~ r^o after the primary special 

collision. 

(ii) The jet (or at least one jet) in the primary event should point in the direction 

of our detector within a A6° uncertainty with cos(A6 |0 ) given by eq.( [148 ) 



Jointly these two requirements will reduce the background by 

3 x 1 -8 — 1 0~ 4 1 1 
(r R o/At H ) (AQ/An) (~ ^-^ ) (± - -) ~ (5 x lO" 10 - 10" 4 ) (150) 

where 

AQ = 2tt[1 - cos(A# )] = 2tt(1 - h ■ r ) = 2n(- -) (151) 
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This background reduction may suffice for facilitating an R° research. 

Finally if the hadronic showers in the special events which do correlate with 
decays in the detector D extend further into the calorimeters than the average 
hadronic (non,b,c) jets, as we would expect on the basis of the R° evolution, we will 
have an additional beautiful confirmation! 

Needless to say all this very qualitative arguments need verification via detailed 
Monte-Carlo studies with g jets and evolving R° particles. 

The whole analysis can be repeated for e + e~ collisions at the Z° resonance. We 
need now to specifically look for 4 jet events: e + e~ — > qqgg which may be relatively 
rare occuring in only few percent of the events. However we gain by having a 
much reduced hadronic background and a larger 5t. [The very large e.m. radiative 
background can be effectively avoided by thick lead shieldings around the detector]. 

In summary we have suggest a novel, and in our view optimized, alternative R° 
search method. Instead of searching for R° decays in fixed target neutral beam where 
the putative R°-s are at a considerable disadvantage, we suggest using colliders and 
jets at 90° where R° could be relatively copious. Also by moderating first the R° 
particles a very wide range of R° lifetimes becomes accessible. 

While many of the basic R° features appear in the papers by G.Farrar we included 
them for completeness. The discussion of sections 4-5 has new elements - and in 
particular the possibility of having longer attenuation path for R° in matter as 
compared with neutrons is discussed at length. 

Negative results in the experiments suggested cannot unequivocally exclude the 
eight gluino hypothesis. Though we argued for R°N cross sections smaller than iV — 
N cross sections the reverse may be true and the R° will then migrate less than our 
estimates making it more difficult to separate R° decays from the "normal" hadronic 
activity. Also for r^o ~ 3 x 10~ 9 — 3 x sec direct searches for decay in flight 
may be more efficient. 

Eventually many different lines of research will converge to yield a definitive 
verdict on the light gluino hypothesis. At present however the R° may exist. It will 
be a great pity if it will be missed due to lack of enthusiasm and possibly for not 
trying the right optimal approach to discovering it. 
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